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Abstract 

Background: Toxin-antitoxin (TA) systems are commonly found in bacteria and Archaea, and it is the most 
common mechanism involved in bacterial programmed cell death or apoptosis. Recently, MazF, the toxin 
component of the toxin-antitoxin module, has been categorized as an endoribonuclease, or it may have a function 
similar to that of a RNA interference enzyme. 

Results: In this paper, with comparative data and phylogenetic analyses, we are able to identify several potential 
MazF-conserved motifs in limited subsets of foodborne pathogens and probiotic strains and further provide a 
molecular basis for the development of engineered/synthetic probiotic strains for the mitigation of foodborne 
illnesses. Our findings also show that some probiotic strains, as fit as many bacterial foodborne pathogens, can be 
genetically categorized into three major groups based on phylogenetic analysis of MazF. In each group, potential 
functional motifs are conserved in phylogenetically distant species, including foodborne pathogens and probiotic 
strains. 

Conclusion: These data provide important knowledge for the identification and computational prediction of 
functional motifs related to programmed cell death. Potential implications of these findings include the use of 
engineered probiotic interventions in food or use of a natural probiotic cocktail with specificity for controlling 
targeted foodborne pathogens. 
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Background In addition, according to the World Health Or- 

Foodborne illnesses continue to be an important public ganization (WHO), antibiotic overuse in food animal 

health concern in developing, as well as in developed production is a major contributor to the emergence of 

countries, thus prevention of foodborne illness outbreaks antibiotic resistant foodborne pathogens [1]. The use of 

through effective and novel interventions should be antibiotics in food animals for growth promotion and 

given priority. The U.S. Public Health Service has identi- treatment disrupts the normal beneficial commensal 

fied ten important foodborne pathogens causing human bacterial microflora in the animal intestinal tract [2-6]. 

illnesses, including pathogenic strains of Escherichia coli, Recently, the human and chicken gut microbiome pro- 

Salrnonella, Listeria, Clostridium botulinum, Shigella, jects [7-12] have shed new light on the existence of a 

and Campylobacter, which are associated with more bacterial phylogenetic core' [13] consisting of a wide di- 

than 250 known foodborne diseases (http://www3.niaid. versity of gastrointestinal bacteria by using new tech- 

nih.gov/topics/foodborne/default.htm). nologies such as next generation sequencing, 16S rRNA 

screens, metagenomics, and metaproteomics. Healthful, 
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Bacillus clausii, Bacillus pumilus, Lactobacillus acidoph- 
ilus, Lactobacillus reuteri, Lactobacillus rhamnosus GG 
(ATCC 53103), Bifidobacterium infantis, Saccharomyces 
boulardii, Lactobacillus ruminis, Lactobacillus johnsonii 
str. NCC 533, and many others. These known probiotic 
strains have been used as dietary supplements, as treat- 
ments for illnesses caused by foodborne pathogens, and 
for disease prevention. Use of probiotic strains not only 
reduce invasion by bacterial pathogens, but also restore 
and maintain an optimal balance of healthy commensal 
bacteria in the human gut via production of antimicro- 
bials [14-23]. 

One of the major mechanisms recognized as being re- 
sponsible for apoptosis, or programmed cell death, and 
production of toxic metabolites in bacteria is through 
the regulation of a wide variety of bacterial toxin- 
antitoxin modules [24-26], such as MazE/MazF, a 
chromosomal toxin-antitoxin module [27-30], plasmid- 
encoded parD [31-33], chpBIK [26,34], relBE [35,36], 
and the PhoQ-PhoP system [37-40]. MazEF is one of the 
most well-studied toxin-antitoxin (TA) systems and has 
been found on the chromosomes of many bacteria. The 
MazF protein has been recently categorized as an endor- 
ibonuclease [41,42] or as a type of RNA interference en- 
zyme [43]. The link between this TA system and 
quorum sensing has also been explored recently through 
a small pentapeptide (NNWNN) called the Extracellular 
Death Factor (EDF) [44]. This small peptide motif (such 
as NNWNN) is known to be an extra-cellular death fac- 
tor in E. coli and other bacterial species. The necessity of 
an "extracellular death factor" (EDF) or "cell death fac- 
tor" (CDF) via MazEF-mediated cell death is a popula- 
tion phenomenon requiring the activation of quorum- 
sensing (cell-to-cell signaling) peptides in bacteria. High 
cell density was found to be associated with elevated 
concentrations of EDF, and the presence of EDF resulted 
in MazF-induced cell death [44]. From a food safety and 
public health perspective, use of EDF or a similar strat- 
egy may be used in place of antibiotics, resulting in less 
usage of antibiotics. We also noticed in one very inter- 
esting study that the induction of toxin MazF and the 
use of antibiotic share a similar mechanism by inhibiting 
the transcription and/or translation of the MazE anti- 
toxin [45]. 

It has been theorized that there is "one toxin for one 
antitoxin" and interestingly MazF, in some bacteria, 
exhibits a selective inhibition of ribosomes and mRNAs 
[43,46]. Numerous strains of probiotic bacteria, such as 
Lactobacillus spp., have been reported to produce anti- 
microbial agents [47], such as bacteriocins, that inhibit 
or kill closely-related species, or even different strains of 
the same species through the inhibition of transcription 
and translation by receptor binding. The antimicrobial 
activities of bacteriocins are due to a heterologous 



subgroup of ribosomally synthesized cationic peptides 
[48]. Nisin, a polycyclic antibacterial peptide 34 amino 
acid residues long, is one of the most studied bacteriocins 
and is produced by many strains of Lactococcus lactis. It 
has been approved by the FDA for use as a food preser- 
vative, and certain probiotic Lactobacillus strains have 
been thoroughly studied and evaluated in vitro and 
in vivo. For example, L. reuteri controls diarrhea in chil- 
dren and suppresses the growth and pathogenicity of 
harmful foodborne pathogens such as Salmonella, E. coli, 
Staphylococcus, and Listeria [49,50] . L. casei GG has been 
used to treat Clostridium difficile infections and to reduce 
intestinal permeability [51-54]. L. reuteri is known to 
produce a broad-spectrum antimicrobial agent, reuterin, 
composed of the natural metabolic compound 3-hydro- 
xypropionaldehyde, which has been used on the surface 
of sausages to inhibit growth of harmful bacteria and fungi 
[15,16]. However, the molecular mechanisms underlying 
the effectiveness of individual probiotic strains have not 
been systematically studied and characterized. Several po- 
tential mechanisms of action, including their ability to 
generate diverse natural toxic metabolites, lactic acid, and 
other organic acids, enzymes, vitamins, and hydrogen per- 
oxide, as well as antimicrobial peptides such as nisin, have 
been well-described [20]. 

The work reported herein explores the experimental 
antimicrobial possibilities and/or procedures for (1) ex- 
pression of an engineered, stress-induced recombinant 
secreted fusion gene encoded by MazF and a small 
extracellular cell death factor (CDF) on the surface or 
extracellular space of recombinant probiotic bacteria or 
(2) for potential application of a cocktail of natural pro- 
biotic strains via experimental in vitro selection to con- 
trol foodborne pathogens for improving the safety and 
quality of foods, as well as improving human health. The 
use of engineered probiotic strains or the natural pro- 
biotic cocktail consisting of mixed probiotic strain popu- 
lations for targeting foodborne pathogens will potentially 
be able to selectively inactivate these pathogens, even in 
complex food matrices. Through gene/motif reshuffling 
[55,56] and computational molecular modeling, this 
engineered and secreted bacterial fusion protein, MazF- 
CDF, which contains an enterokinase (EK) cleavage site 
[57-60] for releasing active CDF and MazF after protein 
secretion, should have a narrow range of applicability, 
limited to inactivating only specific foodborne pathogens 
such as E. coli 0157, Salmonella, Listeria, Campylobac- 
ter, and potentially other human pathogens. This patho- 
gen specificity is due to the fact that the genetically- 
engineered MazF-CDF fusion protein could be modified 
by the inclusion of specific genomic DNA sequences 
from various commensal, as well as pathogenic bacterial 
strains, identified through DNA/protein structure and 
functional comparisons. Additionally, this engineered 
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antimicrobial polypeptide and MazF will be non-toxic to 
beneficial (healthful/probiotic) bacteria, as well as to its 
host that expresses the protein/peptide. Moreover, these 
engineered and secreted CDFs and MazF proteins/pep- 
tides can easily pass through infectious foodborne patho- 
gen cell barriers mediating cell death, and thus could 
potentially reduce the use of deleterious compounds 
such as antibiotics or other harmful chemicals in the 
feed and food industry. 

Results and discussion 

General genetic analysis of probiotics and foodborne 
pathogen genomes 

The genomic information from selected probiotic strains 
and foodborne pathogens is shown in Table 1. It 
revealed little diversity in genomic GC-content in the 
Bifidobacterium genus, while showing an astonishing di- 
versity in the GC-content among Lactobacillus species, 
ranging from 33 to 51.5%. It has been demonstrated that 
genomic GC-content is correlated with a number of fac- 
tors [61], including genome size [62] from species such 
as Lactobacillus, which ranges from 1.8 to 3.4 Mb in 
length. This demonstrates that the GC-content and gen- 
ome size of Lactobacillus genomes may have implica- 
tions related to the biological complexity and adaptation 
of this genus, and could be due to the rate of recombin- 
ation that has been extensively studied in the E. coli gen- 
ome [63]. Knowledge of genetic diversity is fundamental 
to development of synthetic probiotic genomes and/or 
nucleic acid sequence reshuffling strategies. In this 
paper, we demonstrate that the MazF protein is a suit- 
able candidate for the determination of genetic relation- 
ships within sets of MazF proteins in combination with 
with functional motif analysis. 

Ubiquitous existence of MazE/toxin, MazF 

MazEF is a toxin-antitoxin (TA) module widely distribu- 
ted among many bacterial species, including both food- 
borne pathogens and probiotic strains (Table 1, 2), such 
as Lactobacillus acidophilus, Lactobacillus reuteri, Lacto- 
bacillus rhamnosus GG (ATCC 53103), Escherichia coli, 
Selenomonas sputigena, Enterobacter spp., Campylobac- 
ter spp., Citrobacter spp., Thermoanaerobacter tengcon- 
gensis, Pelotomaculum thermopropionicum, Lactobacillus 
casei, Lactobacillus crispatus, Lactobacillus buchneri, 
Bifidobacterium longum subsp. infantis, Clostridium 
botulinum, Clostridium difficile, Vibrio spp., Listeria 
spp., Bacillus spp., Klebsiella variicola, and Salmonella 
enterica. Recent literature and computational analyses 
have shown the presence of many different types of TA 
modules in various localizations, e.g. some TA modules 
have been found within prophages, prophage-like ele- 
ments, and other mobile genetic elements, such as plas- 
mids [31-33]. Due to the existence of varying types of 



toxin-antitoxin modules and possible gene duplication 
events, in this paper we present the above one-to-one 
best matches of chromosomal-encoded mazF orthologs 
and homologs among 75 publically available genome 
sequences of foodborne pathogens and probiotic strains 
(Table 1), and the publically-available databases such as 
NCBI and the Uniprot database (http://www.uniprot. 
org/) in Table 2. 

Phylogenetic analyses and cluster analysis of MazF/ 
antitoxin protein, a growth inhibitor 

The phylogenetic tree in Figure 1 displays the phylogen- 
tic relationships of many well-recognized genera within 
Enterobacteriaceae, including several important food- 
borne pathogens such as pathogenic E. coli, Salmonella, 
Listeria, and Campylobacter, as well as some major pro- 
biotic strains. To build a phylogenetic tree from the data 
in Table 2, the amino acid sequences of all the MazF or 
growth inhibitor proteins were analyzed using the Gen- 
eious software package v5.5.7 with Neighbor-joining 
(NJ) method by applying ClustalW for sequence align- 
ment (Figure 1). Three main clusters (viz., groups 1, 2, 
and 3) are given in Figure 1. At least two representatives 
of potential probiotic strains are listed for each group, 
depending on the complexity of groups. In group 1, the 
potential probiotic, non-pathogenic strains, such as 
Lactobacillus amylovorus, Lactobacillus crispatus, 
Streptococcus, Enterococcus faecium, Enterococcus, Lacto- 
bacillus plantarum, and Lactobacillus rhamnosus, are 
grouped with the foodborne pathogens E. coli, Vibrio 
vulnificus, and Vibrio cholerae. In group 2, the food- 
borne pathogens Enterobacter, Campylobacter upsalien- 
sis, Klebsiella variicola, Salmonella enterica, Shigella 
flexneri, Shigella dysenteriae, and Citrobacter rodentium 
were shown to be genetically closer to Bacillus selenitire- 
ducens, Bacillus halodurans, and Enterococcus faecalis 
In group 3, some probiotic strains, such as the Bacillus, 
Lactobacillus, Leuconostoc, and the Bifidobacterium gen- 
era are categorized along with other major foodborne 
pathogens, such as Clostridium difficile, Listeria monocy- 
togenes, Listeria grayi, and an emerging foodborne 
pathogen Pediococcus pentosaceus. 

Phylogenetic identification of bacterial toxin MazF 
protein motifs and the relationship between gene 
structure and phylogenetic classification 

In Figure 2, it is shown that the number of candidate 
motifs are slightly different in each group, but with a 
high degree of amino acid sequence variability within 
conserved "hot spots" between groups 1, 2, and 3. To 
determine which motifs are the best candidates for the 
engineered MazF construction (discussed in next sec- 
tion) in each individual group, recent studies [43] sug- 
gest that the N-terminal (the first 80 amino acids in 
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Table 1 Genomic information of selected microorganisms considered as potential probiotic strains and some major 
food-borne pathogens used in this study 



Organism 


Strain 


Pathotype 
/or others 


GenBank Accession 
/Assembly 


GC 

(%) 


Genome 
Size 
(~ Mb) 


C~ nrr\n\ilnhnrfor i mcnlioncic 
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H I .D 


1 7 
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I ML. UUJ7 I Z./ 
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1 8 
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I Mfi 1 R31 1 
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39 1 
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LULLUUULIHUb yUbbt.ll 


A I L.L. jjjZj 
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3^ 3 

JJ.J 


1 Q 
I .y 
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pi UUIULIL 
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fin ^ 
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I .y 
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RifiHnhnrforii im nnimnlic e/ /hen Inrfic 
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AD01 1 
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Bifidobacterium dentium 
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2.6 


Bifidobacterium longum subsp. infontis 
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NC_01 1593.1 
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2.8 



Yan et al. Cell & Bioscience 201 2, 2:39 
http://www.cellandbioscience.eom/content/2/1/39 



Page 5 of 1 3 



Table 1 Genomic information of selected microorganisms considered as potential probiotic strains and some major 
food-borne pathogens used in this study (Continued) 
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LULLUUULIHUb [JIUI ILUI UI 1 1 


JDM1 


pi UUIULIL 


Mr 01 7QR4 1 
I ML. U I Z.yO'-r. I 
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Enterobacter sp. 


638 


pathogenic 


NC_009436.1 


52.9 


4.7 


Vibrio vulnificus 


CMCP6 


pathogenic 


NC_004459.2 


46.7 


5.1 


Escherichia coli 


E24377A 


pathogenic 


NC_009801.1 


50.6 


5.2 


Citrobocter rodentium 


ICC168 


pathogenic 


NC_013716.1 


54.6 


5.4 


Klebsiella variicola 


At-22 


pathogenic 


NC_01 3850.1 


57.6 


5.5 


Bacillus megaterium 


QM B1551 


probiotic 


ASM2582v1 


38.0 


5.5 


Escherichia coli 01 57:H7 


Sakai 


pathogenic 


NC_002695.1 


50.5 


5.6 


Bacillus thuringiensis serovar israelensis 


ATCC 35646 


Bio agents 


ASM16769v1 


35.0 


5.9 


Lachnospiraceae bacterium 


3_1_57FAA_CT1 


Non-pathogenic 


Lach_bact_3_1 _57FAA_CT1 _V1 


NA 


7.7 



*: assembly accession. 



Figure 2) of the MazF protein are the most suitable 
motifs. The other motifs in the C-terminal might be re- 
ferred to as 'incorrect' motifs. There are three criteria 
for evaluating the remaining motifs as "incorrect and 
without biological significance", and this is related to: 
the mean hydrophobicity, identity, and the gene struc- 
ture. In Figure 2, in the analysis for the phylogenetic 
identification of bacterial Toxin MazF protein motifs, 
the mean hydrophobicity sequence and identity were 
computed and compared for each sequence; it is inter- 
esting to find that all the compared MazF proteins have 
a higher degree of mean hydrophobicity, identity, and 
more conserved gene structure among several conserved 
amino acid regions, particular in the N-terminal. The se- 
lection of potential functional MazF motifs is discussed 
in the next section. The functional importance of the 



mean hydrophobicity has also been discussed to involve 
mRNA and protein degradation and slower translation 
of mRNA for disordered proteins (http://employees. 
csbsju.edu/hjakubowski/classes/ch331/protstructure/ 
olprotfold.html). 

Bacterial probiotic cocktail strains 

It is important to note that the benefits of probiotics are 
strain-specific [64,65]. A commercial product, VSL#3, 
developed by Sigma-Tau Pharmaceuticals, Inc., provides 
a mixture of probiotic bacteria (B. breve, B. infantis, B. 
longum, L. acidophilus, L. bulgaricus, L. casei, L. plan- 
tarum, and Streptococcus thermophilus) to help protect 
GI tract integrity [66]. In this study, the bacterial pro- 
biotic cocktail strains we propose would be comprised 
of all representatives of groups 1, 2, and 3 shown in 
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Table 2 The genetic characterization of the transcriptional modulator MazF, a chromosomal cell death factor from 
potential probiotic strains and some major food-borne pathogens used in this study 



Organism 


Type of 
organisms 


Strains 


Protein accession 
(GenBank) 


Gene annotation 


E. coli 


Pathogen 
Pathogen 
Antibiotic 
resistance 


EDL933 

Sakai 

B185 


NP. 
NP. 
ZP_ 


_289336.1 
_3 11669.1 
.06660634.1* 


toxin ChpA 


Enterobacter sp. 


pathogen 


638 


ABP60743.1* 


transcriptional modulator of 
MazE/toxin, MazF 


Shigella dysenterioe 


pathogen 


Sd197 


YP_ 


.405833.1* 


toxin ChpB 


Citrobacter rodentium 


pathogen 


ICC168 


YP_ 


.003366767.1* 


Toxin component of the ChpB-ChpS 
toxin-antitoxin system 


Vibrio vulnificus 


pathogen 


CECT4999=R99 


YP_ 


.001393091.1* 


growth inhibitor 


Listeria welshimeri serovar 6b str. 

Listeria monocytogenes 
Listeria seeligeri serovar 1 /2b str. 
Listeria marthii 
Listeria gray! 


pathogen 
pathogen 


SLCC5334 
FSL J 1-208 
SLCC3954 
FSL S4-120 

DSM 20601 


YP_ 
ZP_ 
YP_ 
ZP_ 

ZP_ 


.849071.1 
.05296226.1* 
.003464028.1 
.07870171.1 

.07054243.1* 


PemK family transcriptional regulator 
PemK family transcriptional regulator 
transcriptional regulator, PemK family 
toxin-antitoxin system, antitoxin 
component, MazF family 
MazF family toxin-antitoxin system protein 


Clostridium difficile 


pathogen 


630 

OCD-66c26 
CIP 107932 


YP_ 
ZP_ 
ZP_ 


.001089981.1* 

.05273557.1 

.05323891.1 


putative regulator of cell growth 

ni itativp rpmilatnr nf rpll nrnwth 

|*J*J LO LI V tZT 1 UiaLUI Ul ^tril ^JI^yVVLI 1 

putative regulator of cell growth 


Klebsiella variicola 


pathogen 


At-22 
1_1_55 


YP_ 
ZP_ 


.003438577.1* 
.06548086.1 


transcriptional modulator of MazE/toxin, 
MazF pemK; protein PemK 


Salmonella enterica subsp. 
arizonae serovar 


pathogen 


RSK2980 


YP_ 


.001573441.1* 


hypothetical protein SARI 04525 


Campylobacter showae 


pathogen 


RM3277 


ZP_ 


.05364729.1 


addiction module antitoxin, RelB/DinJ family 


Campylobacter concisus 


pathogen 


13826 


YP_ 


.001466677.1 


addiction module antitoxin 


Campylobacter jejuni subsp. 


pathogen 


IA3902 


ADC28395.1 


prevent-host-death family protein 


Campylobacter upsaliensis 


pathogen 


JV21 


ZP_ 


.07894578.1* 


ChpA/MazF transcriptional modulator 


Shigella flexneri 2a str. 
Shigella flexneri 5 str. 


pathogen 


301 
8401 


NP. 
YP_ 


_709188.1* 
.690766.1 


PemK protein 

growth inhibitor, PemK-like, autoregulated 


Cronobacter sakazakii 


pathogen 


ATCC BAA-894 


YP_ 


.001436394.1 


bifunctional antitoxin/transcriptional 
repressor RelB 


Lactobacillus coleohominis 


probiotic 


101-4-CHN 


ZP_ 


.05553821.1* 


regulatory protein 


Bacillus megaterium 


probiotic 


QM B1551 
DSM 319 


YP_ 
YP_ 


.003560763.1* 
.003595503.1 


endoribonuclease EndoA 


Bacillus thuringiensis 
serovar israelensis 


probiotic 


ATCC 35646 
KBAB4 


ZP 
YP_ 


.00740711.1* 
.001643128.1 


MazF protein 
MazF protein 


Bacillus weihenstephanensis 
Bacillus cereus 




AH1 134 


ZP_ 


.03233130.1 


PemK family 


Coprobacillus sp. 


probiotic 


29_1 


ZP_ 


.08009764.1* 


PemK family transcriptional regulator 


Enterococcus faecal is 
Lactobacillus plantarum 


probiotic 
probiotic 


V583 

TX0102 

TX2134 

HH22 

WCFS1 


NP_81 4592.1 
EFQ1 1602.1* 
ZP_07557645.1* 
ZP_03985009.1 
NP_786238.1* 


PemK family transcriptional regulator 
Toxin-antitoxin system, toxin compoment, MazF 
Toxin-antitoxin system, toxin compoment, MazF 
PemK family transcriptional regulator 
cell growth regulatory protein 


Lactobacillus rhamnosus 


probiotic 


LMS2-1 
Lc 705 


ZP_ 
YP_ 


.04441477.1* 
.003175134.1 


cell growth regulatory protein 
transcriptional modulator of MazE/toxin, MazF 


Lactobacillus johnsonii 


probiotic 


ATCC 33200 


ZP_ 


.04007131.1* 


PemK family growth inhibitor 


Lactobacillus crispatus 


probiotics 


CW-05 


ZP_ 


.07789240.1* 


ppGpp-regulated growth inhibitor 


Lactobacillus gasseri 


probiotic 


ATCC 33323 

224-1 

MV-22 


YP 
ZP_ 
ZP_ 


.815460.1* 

.06262236.1 

.07711551.1 


toxin-antitoxin system, toxin component, 
MazF family 


Lactobacillus easel 


probiotic 


BL23 


YP_ 


.001986080.1 


Cell growth regulatory protein 


Lactobacillus amylovorus 


probiotic 


GRL1118 


YP_ 


.005854914.1* 


transcriptional modulator of MazE/toxin MazF 
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Table 2 The genetic characterization of the transcriptional modulator MazF, a chromosomal cell death factor from 
potential probiotic strains and some major food-borne pathogens used in this study (Continued) 



Vibrio choleroe 


probiotic 


1587 

NCTC 8457 


ZP_01 95061 1.1* 
ZP_0 1969676.1 


transcriptional modulator of MazE/toxin, MazF 


Pediococcus pentosaceus 


probiotic 


ATCC 25745 


YP_805020.1* 


toxin-antitoxin addiction module toxin 
component MazF (an endoRNAse) 


Bacillus amyloliquefaciens 


probiotic 


Y2 


YP_006327269.1* 


Endoribonuclease 


Bacillus pumilus 


probiotic 


SAFR-032 
ATCC 7061 


YP_00 1485694.1* 
ZP_03054587.1 


PemK family growth inhibitor 
endoribonuclease EndoA 


Lactobacillus sake! subsp. sakai 


probiotic 


Sakei 23K 


YP_396224.1 


DNA-binding protein PemK family 


Leuconostoc mesenteroides 
subsp. mesenteroides 
Leuconostoc mesenteroides 
subsp. cremoris 


probiotic 


ATCC 8293 
ATCC 19254 


YP_81 9271.1* 
ZP_0391 3232.1 


toxin-antitoxin addiction module toxin 

component MazF 

PemK family growth inhibitor 


Enterococcus faecal is 


probiotic 


TX0102 
TX0031 


EFQ1 1602.1* 
EFT96737.1 


toxin-antitoxin system, toxin component, 
MazF family 


Bifidobacterium animalis 
subsp. lactis 


probiotic 


HN019 


ZP_02964075.1* 


transcriptional modulator of 
MazE/toxin, MazF 


Bifidobacterium bifidum 


probiotic 


NCIMB41171 


ZP_03645723.1 


RelB antitoxin 


Bifidobacterium longum 
subsp. longum 
Bifidobacterium breve 


probiotic 


JDM301 
DSM 20213 


YP_003660624.1* 
ZP_06596513.1 


RelB antitoxin 

toxin-antitoxin system protein 


Bifidobacterium bifidum 


probiotic 


S17 


YP_003938101.1 


hypothetical protein with RelB 
antitoxin domain 


Lactobacillus easel 


Probiotic 


ATCC 334 


YP_807723.1 


toxin-antitoxin addiction module toxin 
component MazF (an endoRNAse) 


Lactobacillus easel 


Probiotic 


BL23 


YP_00 1988635.1 


Growth inhibitor 


Lactobacillus paracasei subsp. 
paracasei 


Probiotic 


ATCC 25302 


ZP_03963082.1* 


PemK family transcriptional regulator 


Lactobacillus easel str. 


probiotic 


Zhang 


YP_003789562.1 


toxin-antitoxin addiction module toxin 
component MazF 


Lachnospiraceae bacterium 


probiotic 


3_1_57FAA_CT1 


ZP_08609193.1* 


hypothetical protein HMPREF0994_05199 


Leuconostoc citreum 


Probiotic 


KM20 


YP_001 727503.1* 


growth inhibitor 


Lactobacillus reuteri 


probiotic 


DSM 20016 
JCM 1112 
MM2-3 


YP_001 270863.1* 
YP_001 841 242.1 
ZP_03847756.1 


transcriptional modulator of MazE/toxin, MazF 
hypothetical protein LAR_0246 
PemK family growth inhibitor 


Streptococcus sp. oral 
taxon 071 str. 


probiotic 


73H25AP 


ZP_07458552.1* 


ChpA/MazF transcriptional modulator 


Bacillus selenitireducens 


probiotic 


MLS10 


YP_003700699.1* 


transcriptional modulator of MazE/toxin, MazF 


Bacillus halodurans 


probiotic 


0-125 


NP_244588.1* 


ppGpp-regulated growth inhibitor (ChpA/MazF) 


Enterococcus faecium 


probiotic 


DO 

TX0133a04 
TX0133C 


ZP_0571 4797.1* 

ZP_07845579.1 

ZP_07850444.1 


PemK family protein 

toxin component, MazF family 

toxin component, MazF family 



*: data are presented in Figure 1. 



Figure 1. The principal basis behind the composition of 
these probiotic cocktail strains is the assumption that a 
combination of organisms might be more effective than 
the application of a single strain, which potentially 
could suppress many foodborne pathogens, such as the 
E. coli and vibros in Group 1, several foodborne patho- 
gens, such as Enterobacter, Klebsiella variicola, Salmon- 
ella enterica, Shigella flexneri, Shigella dysenteriae, 
Citrobacter rodentium, and Campylobacter upsaliensis, 
(one of the most common Campylobacter strain found in 
people with diarrhea in Group 2), Clostridium difficile, 



Listeria monocytogenes, and Listeria grayi, in Group 3 of 
Figure 1. Table 1 shows an astonishing diversity in the 
genomic GC-content and genome size among Lactobacil- 
lus species and their diverse distribution in all groups 
within Figure 1, which indicates a potential to further 
identify other closely-related Lactobacillus species (not 
listed in Table 1) into the three previously-described 
groups. The hypothesis underlying our approach is that 
probiotic strains found within the same group with food- 
borne pathogens will have a reasonable degree of genetic 
and molecular phylogenetic compatibility and could 
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1. Lactobacillus amylovorus | YP_005854914.1 

2. Lactobacillus crispatus |ZP_07789240.1 

3. E. coli |ZP_06660634.1 

4. Vibrio vulnificus |YP_001393091.1 

5. Vibrio cholerae |ZP_01950611.1 

6. Streptococcus |ZP_07458552.1 

7. Enterococcus faecium |ZP_05714797.1 

8. Enterococcus |EFQ11602.1 

9. Lactobacillus plantarum |NP_786238.1 

10. Lactobacillus rhamnosus | ZP_04441477. 1 - 

11. Enterobacter | ABP60743.1 

12. Campylobacter upsaliensis |ZP_07894578.1 

13. Bacillus selenitireducens | YP_003700699.1 

14. Bacillus halodurans |NP_244588.1 

15. Enterococcus faecalis |ZP_07557645.1 

16. Klebsiella variicola |YP_003438577.1 

17. Salmonella enterica |YP_001573441.1 

18. Shigella flexneri |NP_709 188.1 

19. Shigella dysenteriae | YP_405833.1 

20. Citrobacter rodentium | YP_003366767. 1 - 

21. Bifidobacterium lactis | ZP_02964075. 1 

22. Lactobacillus johnsonii |ZP_04007131.1 

23. Lactobacillus gasseri |YP_815460.1 

24. Lachnospiraceae bacterium |ZP_08609193.1 

25. Lactobacillus coleohominis |ZP_05553821.1 

26. Coprobacillus |ZP_08009764.1 

27. Clostridium difficile |YP_001089981.1 

28. Listeria monocytogenes |ZP_05296226.1 

29. Listeria grayi |ZP_07054243.1 

30. Bacillus thuringiensis |ZP_00740711.1 

31. Bacillus megaterium | YP_003560763.1 

32. Bacillus amyloliquefaciens | YP_006327269.1 

33. Bacillus pumilus | YP_001485694.1 

34. Lactobacillus sakei |YP_396224.1 

35. Lactobacillus paracasei |ZP_03963082.1 

36. Pediococcus pentosaceus |YP_805020.1 

37. Lactobacillus reuteri | YP_001270863.1 

38. Leuconostoc mesenteroides |YP_819271.1 

39. Leuconostoc citreum | YP_001727503.1 

Figure 1 Phylogenetic analyses of Enterobacteriaceae based on MazF (toxin), a growth inhibitor. The NCBI accession numbers are shown 
after each taxon name. 
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bridge a relationship similar to a "symbiosis" of entities, 
including exchanging toxin/antitoxin molecules among 
the probiotic and pathogenic strains. Lactobacillus spe- 
cies are known to produce antimicrobial substances, in- 
cluding bacteriocins, lactic acid, and hydrogen peroxide. 
The MazEF toxin component may provide a basis for 
bacterial growth inhibition within the same group (Fig- 
ure 1 and 2). Therefore, this toxin-antitoxin module may 
have great potential to inhibit the growth of potentially- 
pathogenic bacteria through a possible competitive ex- 
clusion due to selective inhibition [46]. Figures 1 and 2 
list all possible cocktails of these probiotic strains. In 
reality, a foodborne outbreak is more likely to be asso- 
ciated with one particular foodborne pathogen in par- 
ticular foods. For example, there is a low incidence of 
Campylobacter in ground beef and pork, while Cam- 
pylobacter is the major foodborne pathogens associated 
with poultry, therefore, a single food-borne pathogen 
with the application of mixed probiotics will be consid- 
ered initially. 



Molecular recombination techniques: construction of 
genetically engineered synthetic probiotic strains 

Figure 3 details an engineered probiotic strain bearing a 
recombinant plasmid containing a stress-induced pro- 
moter, followed by an in-frame gene fusion accom- 
plished by fusing an appropriate signal peptide, a 
functional cell death peptide/factor (CDF), an enteroki- 
nase (EK) binding site, and a genetically-modified engi- 
neered MazF gene, which will be constructed and 
transformed into the probiotic bacteria. In this recom- 
binant probiotic strain, environmental stress, relevant to 
food environments will trigger gene expression of this 
fusion protein under an environmental stress-inducible 
promoter. The signal peptide directs the encoded fusion 
protein to the extracellular space of the engineered pro- 
biotic strains. The signal peptide depleted fusion pro- 
teins will be cleaved by a biliary tract enterokinase 
directly into the digestive system to release the 
functionally- active CDF and MazF. This event will occur 
based on the ability of probiotic strains to form stable, 
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Figure 2 ClustalW output (default settings) for the aligned amino acid sequences of MazF. The pdf images of the alignments were 
generated using Geneious v5.5.7. Conserved amino acids are highlighted in colors; the numbers (1 to 39) in this Figure correspond to the taxon 
names in Figure 1. 



non-infectious biofilm-type aggregates, which may attach 
with other bacteria, including foodborne pathogens, in 
the urogenital and intestinal tracts [22]. These active 
species-specific CDF and MazF proteins will bind and 
pass through "unfriendly" bacterial cell surface barriers 
into the cytoplasm. These processed CDFs and MazF 



proteins will selectively inactivate and/or inhibit proteins 
involved in cell survival and induce the synthesis of more 
cell death-related proteins with activity against food- 
borne pathogens, eventually controlling, inhibiting, or in- 
activating "unfriendly" cells. Although antitoxin MazE 
could reverse the bacteriocidal effect of the overexpressed 



"Friendly" Probiotic Bacteria 



L. acidophilus, L salivarius, 
L. easel, L. thermophilus, 
B. hifulum, B. longum, ete 



•Unfriendly" Food-borne Pathogens 



Plasma mmbiane 1«<«"™' W«™««'<»> 



Enterokinase (EK) £ 
from biliary tract EK CDF 



Salmonella, Listeria, 
E.coli 0157, 
Campylobacter, etc 




i— i CDF and mazF wi 1 1 
selectively inhibit or 
inactivate genes 
„ involved in cell 
survival and cell death 
related proteins 



i 



J > Cell death 



Figure 3 A Schematic representation of the engineered, secreted CDF and MazF-mediated cell death during food processing and 
storage. SP: promoter, CDF: cell death factor/peptide, EKB: enterokinase binding site. 
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MazF, MazE cannot impede the downstream cascade 
already initiated by MazF at early stages of the MazF- 
mediated cascade [42]. 

The overall hypothesis for this experiment is graphic- 
ally presented in Figure 3. The engineered MazF gene se- 
quence will be designed based on the genomic sequences 
of all publically- available foodborne pathogenic strains by 
using reasoned random gene biosynthesis and/or genuine 
gene reshuffling to rapidly combine functions and prop- 
erties of parental genes for the development of improved 
gene specificity and generality, molecular modeling (tran- 
scription factor binding site identification, etc.), and sys- 
tems biology technologies/tools. 

Conclusions 

Survival of foodborne pathogens in cultures or in animal 
GI tracts may be very genus- or species-specific. Data 
presented in this paper can be explored to develop ef- 
fective intervention strategies applied directly during 
food processing and preparation, as well as in the animal 
feed supply, which may lead to an overall reduction in 
use of antibiotic growth promoters (AGP) throughout 
the world. Recent research in molecular biology and 
genomics has provided potential applications of pro- 
biotic strains as dietary supplements, which could re- 
place AGP in animal diets or as biotherapeutic agents in 
cases of antibiotic-associated diarrhea in travelers and in 
childhood diarrhea and other bacterial gastrointestinal 
illnesses. Experiments relating to this potential probiotic 
application may reveal a further greater range of poten- 
tial benefits. For many of these potential benefits, cur- 
rent research is limited, and only preliminary results are 
available. All effects can only be attributed to the indi- 
vidual strain(s) tested. Testing of a specific supplement 
may not be extrapolated to benefits from any other 
strain of the same species, and testing results do not 
imply that comparable benefits will be imparted from 
other LAB (or other probiotics). In this study, we have 
computationally explored several potential intervention 
strategies to control foodborne pathogens, either by 
using a cocktail of probiotic strains or an engineered 
probiotic strain. 

The inhibition of pathogens by probiotic strains is 
mainly due to the production of antibacterial peptides 
[67], the release of short-chain fatty acids, or reduction 
of the pH within the lumen [68,69] by the production of 
organic acids or by decreasing pathogen adherence to 
intestinal epithelial cells [70]. Therefore, the benefits of 
probiotics could be very strain- or species-specific, and 
probiotic strains may rely on different mechanisms to 
suppress growth, attachment, or other metabolic pro- 
cesses, inherent to pathogenic bacteria. Moreover, the 
optimal effects of probiotic strains may involve the sim- 
ultaneous use of more than one strain. Our contention 



in this paper is not experimental proof but, rather, a 
clear scientifically-backed hypothesis in the form of a 
detailed accompanying method that multi-probiotic 
strain composites with diverse genetic backgrounds may 
complement [71] one another as vectors of competitive 
exclusion and, therefore, could maximize the potential 
to inhibit an array of common foodborne pathogens [72] 
in the gastrointestinal tract of humans or livestock, as 
well as in foods and animal feed. 

The use of probiotic bacteria with the ability to pro- 
duce CDFs and engineered MazF to selectively inactivate 
pathogens is a novel approach to controlling pathogens 
in foods, and possibly treating human infections. A 
number of studies suggest that this project could have 
practical significance and be a potentially new approach 
for the development of novel and cost-effective food 
safety intervention technologies for the control of food- 
borne pathogens and improving public health [73-75]. 

The approaches described above represent a first at- 
tempt to describe a systematic approach or method to 
test the hypothesis that "friendly" bacteria can be used 
to inactivate or inhibit pathogens in food based on ex- 
pression of MazF. There are many specific cell death fac- 
tors that may be associated with bacterial programmed 
cell death and multi-cellular behavior mechanisms in 
foodborne pathogens. Through computational modeling, 
remodeling, genetic recombination, or further gene re- 
shuffling, and exploring experimental approaches, it may 
be possible to evaluate and elucidate more effective 
CDFs and MazF to be used for controlling foodborne 
pathogens, which will ultimately result in a reduction in 
the use of antimicrobial compounds in humans and ani- 
mals, as well as during food processing and storage. 

Materials and methods 

Genomic sequences 

All of the genome and gene sequences examined in this 
study (Table 1 & 2) are available in GenBank (http://www. 
ncbi.nlm.nih.gov/genbank/GenbankOverview.html). 

Identification and analysis of bacterial toxin-antitoxin 
modules, MazE/MazF 

The prediction accuracy of the best chromosomal- 
encoded MazF orthologs among relatively distinct gen- 
ome strains is critical for the performance of molecular 
phylogenetic analysis. We used a sequential BLAST 
workflow based on pairwise comparison by applying 
either an E. coli programmed cell death toxin MazF 
(Genbank accession: ZP_06660634.1), an endoribonu- 
clease MazF, a MazF protein from Vibrio cholerae 
(ZP_01950611.1), or a PemK family transcriptional 
regulator protein (ZP_05296226.1) from Listeria 
monocytogenes to perform a BLASTP homology search 
of these combined protein sequences with all 75 
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publically available foodborne pathogen and probiotic 
strain genomic sequences (presented in Table 1) publi- 
cally available databases such as NCBI and Uniprot data- 
base (http://www.uniprot.org/). BLAST search results 
were used to list MazF or MazF-like proteins from 75 dif- 
ferent strains being represented as 39 different species 
(Table 2). 

Phylogenetic analysis and computational identification of 
phylogenetic motifs of the MazF protein 

In order to identify potential functional motifs of MazF 
proteins, phylogentic analyses of MazF proteins were 
conducted by applying the embedded multiple sequence 
alignment ClustalW program in the Geneious software 
package v 5.5.7 [76,77] with the Neighbor Joining method. 
ClustalW output for the aligned amino acid sequences 
and the pdf images of the alignments were generated 
using the Geneious software package v 5.5.7. 
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